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Abstract 

A new uniform asymptotic approximation for the Wigner 6j sym- 
bol is given in terms of Wigner rotation matrices (<i-matrices). The 
approximation is uniform in the sense that it applies for all values of 
the quantum numbers, even those near caustics. The derivation of 
the new approximation is not given, but the geometrical ideas sup- 
porting it are discussed and numerical tests are presented, including 
comparisons with the exact 6j-symbol and with the Ponzano-Regge 
approximation . 

1 Introduction 

The Wigner 6j-symbol is used in the recoupling of three angular momenta, 
and finds many applications in atomic, molecular and nuclear physics. These 
are explained in standard references on angular momentum theory [Tj, El IBIS]. 
For example, the 3nj-symbols and their asymptotic properties are central to 
certain algorithms for the calculation of scattering amplitudes in three-body 
interactions (De Fazio et al[5J, Anderson and Aquilanti[6j). These methods 
make use of the relationship between the 3nj-symbols and discrete orthogonal 
polynomials (Aquilanti et al[7J El [9] and references therein). 
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The 6j-symbol possesses a remarkable semiclassical approximation, first 
obtained by Ponzano and ReggefTU] through some inspired guesswork, that 
is linked in a highly symmetrical manner to the geometry of a tetrahedron in 
three-dimensional space. This formula was first proven rigorously by Schulten 
and Gordon[llJ, who also gave practical means of computing the 6j-symbol 
by recursion relations. More recently the 6j-symbol has attracted attention 
for the role it plays in quantum gravity, which has led to more geometri- 
cal treatments of its asymptotic properties. References in this area include 
Roberts [T2] and Charles [T3]. 

The formula of Ponzano and Regge has the usual properties of primitive 
semiclassical approximations, for example, it diverges at the classical turning 
points (the caustics). Since the 6j-symbol is only defined for discrete values 
of the quantum numbers, it is unlikely to fall exactly on a caustic, but it can 
come very close, and the Ponzano-Regge formula does not provide a good 
approximation for such values. Thus there is interest in uniform approxima- 
tions that do not suffer from caustic singularities and that are valid over as 
wide a range of quantum numbers as possible. 

In addition to proving the Ponzano-Regge formula for the 6j-symbol, 
Schulten and Gordon also provided uniform asymptotic approximations of 
the Airy-function type that are valid in a region passing through a turning 
point. The 6j-symbol, however, has two turning points when one of the j's 
is varied and the others held fixed, in a manner reminiscent of an ordinary 
oscillator in one dimension. The Airy-function uniform approximation can 
only cover one of these at a time, and if two Airy-function approximations are 
used, then they do not match smoothly in the middle. This suggests that a 
uniform approximation of the Weber-function type (one based on harmonic 
oscillator wave functions) should be used, that would cover both turning 
points at once. 

It turns out, however, that this cannot be done, in general. An idea of 
why this is so is given by Fig. [TJ which shows the 6j-symbol and the Ponzano- 
Regge approximation as a function of ji2 for certain values of the other five j's 
(notation is explained by ([T]) below). The sticks in the diagram give the exact 
values of the 6j-symbol, while the curve is the Ponzano-Regge approximation. 
Although the 6j-symbol is only defined for discrete values of j'12, the Ponzano- 
Regge approximation is defined for continuous values of the parameters, and 
allows us to talk of the phase of the 6j-symbol in an unambiguous manner. 
The Ponzano-Regge approximation shows the expected divergences at the 
caustics or turning points, and the Airy function behavior at the right turning 
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point. In the classically allowed region between the turning points the wave 
length decreases as we move to the left, corresponding to an increase in the 
variable conjugate to j i2 (the angle <pi2, defined in Sec. !2.4b elow). At the left 
turning point the Airy function is not so obvious, because it is multiplied by 
a rapidly oscillating cosine term. This is because the angle 0i2 conjugate to 
ji2 reaches the value ±7r at the right turning point. 

The behavior illustrated in Fig. [JJ cannot be matched by a Weber function 
(harmonic oscillator eigenfunction) , for which both turning points have the 
behavior of an unmodulated Airy function. Related to this is the fact that 
the difference between the Ponzano-Regge phases at the two turning points 
is not of the form (n + |)7r for the parameters shown in Fig. [TJ For har- 
monic oscillator eigenfunctions the difference in the action between the two 
turning points is always of the form (n + ^)tt, a condition that is equivalent 
to the single-valuedness of the WKB wave function (that is, it is the Bohr- 
Sommerfeld quantization condition). One can say that the 6j-symbol has 
a nonstandard matching condition of the two WKB branches at the lower 
turning point for the parameters in Fig. [TJ 

Ultimately there are topological reasons for the failure of the Weber func- 
tion as a standard form for a uniform approximation in cases such as that 
illustrated in Fig. [TJ Uniform approximations are based on a smooth, area 
preserving map between the phase space of the given problem and the stan- 
dard problem. But the phase space of the 6j-symbol is a sphere (we call 
it the 6j'-sphere), and the phase space of the harmonic oscillator is a plane. 
These two spaces cannot be continuously mapped into one another. 

We have noticed, however, that the phase space that arises in the semi- 
classical analysis of the Wigner ci-matrices (rotation matrices) is also a sphere. 
We call this the ci-sphere, to distinguish it from the 6j-sphere. Moreover, the 
classical orbits on the two spheres are topological circles in both cases, and in 
both cases intersecting orbits always intersect in two points, unless they are 
tangent, in which case there is one intersection point. There are also cases in 
which the orbits do not intersect at all, corresponding to classically forbidden 
regions. As is well known, the classical orbits and their intersections provide 
the geometrical framework for the construction of semiclassical (asymptotic) 
approximations. These topological features are the same for both the Gj- 
sphere and the (i-sphere, suggesting that one can be mapped into the other 
by a smooth transformation that takes a pair of orbits on one sphere into the 
pair of orbits on the other. 

We have worked out the details of this mapping and the corresponding 
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uniform approximation. The resulting approximation is smooth and uniform 
over the entire range of quantum numbers ji 2 and j 2 3 (for fixed ji, j 2 , jz and 
ji). For most values of the j's it is more accurate than the Ponzano-Regge 
approximation, certainly near the caustics but also at most other places. We 
have found no places where it is dramatically worse than the Ponzano-Regge 
approximation. 

In this paper we shall present the new uniform approximation itself, as 
well as some of the geometric rationale behind it, which helps considerably in 
understanding the formula and the various regions that it covers. In addition, 
we shall present some numerical tests of the new formula and comparisons 
with the exact 6j-symbol and the Ponzano-Regge approximation. We shall 
not, however, present the details of the derivation. 

The outline of this paper is as follows. In Sec. [2] we present a collection of 
facts about the 6j-symbol and its asymptotic or semiclassical approximation 
(the Ponzano-Regge formula), emphasizing the spherical phase space and its 
geometrical ramifications. In Sec.[3]we present a geometrical treatment of the 
asymptotic properties of the rotation matrices (d-matrices) that emphasizes 
the similarities with the 6j-symbol. In Sec. H]we outline the ideas behind the 
uniform approximation of the 6j-symbol in terms of <i-matrices, we present 
the actual formulas, and we present some numerical tests. Finally, in Sec. [5] 
we make some conclusions. 

2 The 6j/-symbol 

2.1 Quantum Mechanics of the 6j-symbol 

We set ft = 1, so all angular momenta are dimensionless. We label the j's in 
the 6j symbol by 



which is how it would used when recoupling three angular momenta. 

The quantum number ji, i — 1, 2, 3, 4, 12, 23, just gives the magnitude of 
the angular momentum, and does not specify the sign of the operator. For 
example, instead of coupling three angular momenta to obtain a fourth, that 
is, setting J 4 = Ji + J 2 + J3, we can couple four angular momenta with a 
sum of zero, 




(1) 



Ji + J 2 + J 3 + J 4 = 



(2) 
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(effectively changing the sign of J4). This is how we shall regard the recou- 
pling problem in this paper. Usually we will think of ji, i = 1, 2, 3, 4 as given, 
while ji2 and j'23 are variable intermediate angular momenta that result from 
the coupling of the first four. They are the quantum numbers of the squares 
of the operators 

Jl2 = Jl+J2, J23 = J2 + J3- (3) 

With this interpretation, the 6j-symbol in the form ([T]) is proportional 
to the unitary matrix element (J12 1 J23) that takes one from the eigenbasis 
of one intermediate angular momentum (ju) to the eigenbasis of the other 
0/23 )• These bases span the subspace of the product space of four angular 
momenta in which ([5]) holds as an operator equation. We shall denote this 
subspace by Z. According to (|2J), the total angular momentum vanishes on 
Z. The orthonormality relations satisfied by the 6j-symbol (see, for example, 
Edmondsjl] Eq. (6.2.9)) are essentially a statement of the unitarity of the 
matrix <ji 2 |j23>- 

To be defined the 6j-symbol (pQ) must satisfy four triangle inequalities, in 
(ji, 32, J12), U2J3J23), (J3,J4, J12), and (j u j 4 , j 23 ), for example, ju must lie 
between the bounds 

\ji ~ J2I < J 12 < Ji +32, (4) 

in integer steps. For given ji, i = 1,2, 3, 4, these imply that and J23 vary 
between the limits 



Jl2,i 
J23,] 



< 



Jl2 ^ Jl2,max) 
< J23 < J23,max, 



in integer steps, where 



Jl2,min 
J23,min 
Jl2,max 
J23,max 



max(|ji -j 2 \, \ja-j4), 
max(|j 2 \ji ~ 

min(ji + J2,h+k), 
min(j 2 +h,h+k)- 



(5) 



(6) 



The number of allowed j'12 or J23 values is the same, and it is the dimension 
D of the subspace Z as well as the size of the matrix (j'12 1 J23) , 



D — dim Z — jl2,max — Jl2,min + 1 — J23,max ~ J23,min + 1- 



(7) 
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2.2 Classical and Semiclassical Mechanics of the 6j- 
symbol 

The basic reference on the semiclassical mechanics of the 6j-symbol is Pon- 
zano and Regge[10j. We add to their discussion an appreciation of the Gram 
matrix (see Appendix |A|) and the recent realization that the phase space of 
the 6j-symbol is a sphere [13]. 

We shall reserve lower case for quantum numbers as in ([1]), and for 
semiclassical purposes we shall set 

J i=3i + l ( 8 ) 

(with capital J's), for i = 1,2,3,4,12,23. The quantity J, t is interpreted 
as the length of the classical angular momentum vector associated with the 
quantum number j%- The correction 1/2 is a Maslov index, and the manner in 
which it arises in the semiclassical theory of angular momentum is explained 
in Aquilanti et al[14j. 

We shall write J, (in bold face) either for the vector of angular momentum 
operators in quantum mechanics, or for the classical angular momentum 
vector in a classical model. The distinction will be established by context. 
For example, in the classical context, (j2J) is interpreted as a constraint on 
the four classical vectors Jj, i = 1,2,3,4, while is interpreted as the 
definitions of two more classical angular momenta J12 and J23. 

The condition (jSj) is a quantization condition of the Bohr-Sommerfeld 
type, restricting the the classical quantity Jj to discrete values. Some authors 
have viewed (|5J) as an approximation to + l)] 1 ^ 2 , valid when j is large, 
but in fact it is represents the exact eigenvalues of a certain operator for all 
values of j, even j = 0. When properly understood [2], (jSJ) is equivalent to 
the fact that the eigenvalues of the operator J 2 are j(j + 1). 

In purely classical mechanics, however, there is no quantization, and all 
variables take on continuous values. We must allow this in order to view the 
classical phase space. To visualize the phase space of the 6j-symbol, we will 
assume that Jj, i = 1, 2, 3, 4 have any fixed positive values, while J12 and J23 
are variables, the lengths of the vectors J i2 and J23 defined by (jSJ). Then J 12 
and J23 are restricted by classical versions of the triangle inequalities, 

</l2,min < Jl2 ^ ^12,max, 

^23,min — ^23 — ^23,max ; (9) 
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where Ji 2 and J 23 vary continuously between the bounds indicated, and where 
the bounds themselves are given by 

</i2,min = max(|Ji - J 2 |, |J 3 - J 4 |), 

^23,min = max(| J 2 - J3I , | Jl ~ M) , 

^i2,ma X = min(Ji + J 2 , J 3 + J 4 ), 

^23,max = min( J 2 + J3, Jl + J 4) ■ (10) 

These are the bounds of Ji 2 and J 23 on the 6j-sphere. 

Here are two useful theorems. First, if J 23 , m in = J\ — -h or J 2 — J 3 , then 
</i2,max = J3 + Ji, otherwise Ji 2 , max = Ji + J 2 . Second, if Ji 2 , mm = Ji - J 2 or 
J 4 - J 3 , then J 23 ,max = ^2 + ^3, otherwise J 23 ,max = Ji + ^4- 

If the quantization conditions ([S]) hold for J i5 z = 1,2,3,4, then the 
bounds on the continuous variables Ji 2 and J 23 can be expressed in terms of 
the bounds on the quantum number ji 2 and j 23 by 

>-^12,min Jl2,min; >-^12,max Jl2,max 1; 

^23 ,min — J23,miii) ,max — J23,max 

+1. (11) 

Combined with (17j) , these imply 

-D = t/l2,max — ^12,min = ^23,max — </23,min- (12) 

If Ji > 0, i — 1,2,3,4 and if J i2 and J 23 satisfy the triangle inequalities 
IQ, then it is possible to find six classical vectors Jj, i = 1, 2, 3, 4, 12, 23 that 
can be placed end-to-end in subsets of three to create four triangles. The 
triangles are defined by d3J), plus Ji 2 + J 3 + J 4 = and Ji + J 23 + J4 = 0. In 
particular, this can always be done when the six Jj satisfy the quantization 
conditions (jBD for values of ji that are valid in a 6j-symbol, in which case 
Ji > for all six i (because ji > 0) and the areas of the triangles are positive 
(because the ji satisfy the triangle inequalities). 

For some values of the six Jj the four triangles can be fitted together to 
form the four faces of a tetrahedron. This is the classically allowed region 
of the 6j-symbol. If they can, then the signed volume of the tetrahedron is 
given by 

V = -U • (J 2 x J 3 ) = ■ (A 2 x A 3 ), (13) 
o 

where the vectors Aj are defined in fl86l) . The volume is related to the 
nonnegative definite Gram matrix G defined in ( 1ST)) by 

3QV 2 = detG. (14) 
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The tetrahedron is illustrated in Fig. [2J which shows our convention for 
labeling the edges by classical angular momentum vectors. The vectors in 
the figure may be seen to satisfy (j2J) and ([3]). The tetrahedron in Fig. [2] has 
positive volume (V > 0) according to the definition f[T3"[) . 

For other values of the Ji a real tetrahedron does not exist, but a complex 
tetrahedron exists whose edges are complex vectors Jj. These may be chosen 
so that the x- and y-components are purely real and the z-component is 
purely imaginary. This is the classically forbidden region of the 6j'-symbol. 
In this case G is still given by (IHTj) and it is still a real, symmetric matrix, 
but it has one negative eigenvalue. Also, V 2 < and V is purely imaginary. 
Equation (I13p is still valid, however, in terms of the complex vectors Jj. The 
dot product of complex vectors U and V is defined by J2i UiVi (not J2i U*Vi), 
and the length by Jf = Jj ■ Jj (not J* • J,). 

The classically allowed and forbidden regions are illustrated in Fig. [3j 
This figure shows the square region of the Ji 2 -J 23 plane bounded by the 
classical limits ([TO]) for certain fixed values of jj, i = 1,2,3,4. The small 
spots inside the square are the quantized values of J i2 and J 23 , from which 
the quantum numbers ji 2 and j 23 can be extracted by (jSJ). Notice that the 
quantized values of Ji 2 and J 23 always stay at least one half unit away from 
the bounds (fTUl) . The oval curve is the curve V = 0, separating the classically 
allowed from the classically forbidden regions; it is the caustic curve. The 
classically allowed and forbidden regions lie inside and outside the caustic 
curve, respectively. Actually, there are four classically forbidden regions, 
labeled ABCD in the figure. Similar diagrams describe radiative transitions 
in hydrogen (see Fig. 27 of Ref. [T5]). 

The caustic curve consists of real tetrahedra that are flat; it touches the 
square at four points, labeled XYZW in Fig. [3J The geometrical meaning 
of these points and the behavior of the flat tetrahedron as we move around 
the caustic curve are illustrated in Fig. HI At point Y, vectors Ji and J 2 are 
antiparallel, giving Ji 2 its minimum length. As J 2 rotates in a counterclock- 
wise direction, at first Ji 2 grows and J 23 shrinks, as illustrated in the diagram 
Y — > X in the figure. This is a point between Y and X on the caustic line of 
Fig. [31 When J 2 rotates to the angle that causes J 2 and J 3 to be antiparallel, 
then J 23 is at its minimum length and we are at point X. As J 2 continues to 
rotate, J 23 starts to grow again while Ji 2 continues to grow, as illustrated in 
the diagram X — > Z in the figure. This is a point between X and Z on the 
caustic curve. In this manner we may continue around the caustic curve. 
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2.3 The Ponzano-Regge Formula 

Suppose we are in the classically allowed region so a real tetrahedron exists. 
Let ipi, i = 1,2, 3, 4, 12, 23 be the exterior dihedral angles of the tetrahedron 
associated with edge Jj, that is, ipi is the angle between the outward pointing 
normals of the two faces that meet in edge i, so that each ipi lies in the interval 
[0, 7r]. Then the Ponzano-Regge phase is defined by 

$ PR = H^, (15) 

i 

where the sum runs over all six edges, and the Ponzano-Regge approximation 
is 

where the amplitude in the classically allowed region is given by 

^4pr = J- (17) 

'l2ir\V\ 



Given the six quantum numbers ji, an algorithm for determining the 
dihedral angles ipi is the following. The rules we give are equivalent to those 
of Ponzano and Regge and Schulten and Gordon, but stated in terms of the 
diagonalization of the Gram matrix. First we define the six Jj by (j5J and 
then we set up the Gram matrix using (1551) and diagonalize it. If all the 
eigenvalues are positive (or if det G = 36V 2 > 0), we are in the classically 
allowed region and we may proceed. Then we construct the six vectors Jj 
as explained in Appendix |X] we compute the outward pointing normals by 
taking cross products of the vectors spanning the four faces, and finally we 
compute cos^j as the dot products of the outward pointing normals. This 
determines ipi uniquely as an angle in [0, 7r]. This is not the most efficient 
algorithm from a numerical standpoint, since to determine the ipi only the 
dot products of the vectors are needed and not the vectors themselves, but it 
is conceptually clean and has the benefit of allowing one to draw or visualize 
the tetrahedron itself. 

The Ponzano-Regge phase <£>pr is continuous inside the classically allowed 
region, as are the dihedral angles fy. On the caustic boundary all tetrahedra 
are flat so all dihedral angles are either or tt. These angles are continuous 
(hence constant at or tt) on the caustic line between points XYZW, but 
some angles jump discontinuously between and tt at those points. 
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Table 1: The dihedral angles ipi on the segments of the caustic curve bounding 
classically forbidden regions ABCD. There are two possibilities for segments 
B and D, and one for segments A and C. 

The possible values of the ipi on the caustic curve are summarized in Ta- 
ble [TJ The four segments of the caustic curve are identified by the classically 
forbidden region (ABCD) to which they are adjacent. In segments A, B and 
D there are two possibilities, while in segment C there is only one. In seg- 
ment A, the first column applies if Ji2,max = -h + J4 and the second column 
otherwise; in segment B, the first column applies if J23,min = J2 — J3 or J4 — J1 
and the second column otherwise; and in segment D, the first column applies 
if ^23,max = -h + J3, and the second column otherwise. Examples of these 
rules may be seen in Fig. HI The dihedral angles in a flat tetrahedron such 
as the ones labeled Y — > X and X — > Z are for interior segments and tt for 
segments bounding the outside of the plane figure. 

Although some angles tpi are discontinuous at points XYZW, the Pon- 
zano-Regge phase $pr is continuous everywhere on the caustic boundary 
(hence everywhere inside and on the caustic boundary). 

The angles ipi that are tt on a segment of the caustic curve correspond 
to the vectors Jj that lie on the outside of the plane figure, as seen in the 
examples in Fig. @J That is, they correspond to a set of vectors Jj that sum 
to zero. But this implies that the sum of the corresponding ji values is an 
integer, 

% = Ji = inte S er > ( 18 ) 

i 

where the prime means to sum only over i such that ipi on a segment of the 
caustic curve has the value tt. For example, from the first column for segment 
B in Table [1] we have j\ + j'3 + j'12 + J23 = integer. 

In the classically forbidden regions the method of Appendix |X] yields 
vectors Jj whose (x,y,z) components can be labeled as (r,r,i), where r 
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means real and i means imaginary. These vectors have real lengths Jj that 
are fixed by (|S]) and the values of the jj. The cross products Jj x Jj have 
the form (i, i, r) and also have real lengths, which are twice the real areas 
of the faces. Dividing by these we obtain complex unit normals to the faces 
of the form (i,i,r), whose dot products, the cosines of the angles ipi, are 
real. These cosines lie outside the range [— 1,+1], however, indicating that 
the ipi are complex. Since the complex inverse cosine function has multiple 
branches, we must determine the branch. 

A first requirement is that branch chosen for ipi should agree with the 
value of ipi as we approach caustic curve, where ipi is either or tt, depending 
on % and the region ABCD, as shown in Table [TJ If ipi — 0, cosipi = +1 on 
the caustic curve, then cos ipi is real and > 1 in the classically forbidden 
region. In this case we choose ipi = iipi, where ipi = cosh - (cosipi) is real and 
positive. If ipi = ii, cosipi = —1 on the caustic curve, then cosipi is real and 
< — 1 in the classically forbidden region. In this case we choose ipi = n + iipi, 
where ipi = — cosh _1 (— cosipi) is real and negative. We can summarize these 
two cases by 

ipi = sign (cos ?/>j) cosh -1 (| cos^D- (19) 

In spite of the sign and absolute value functions, ipi is a smooth function of 
position in any of the four classically forbidden regions. 

This procedure allows us to determine which classically forbidden region 
(ABCD) we are in, for if we note the signs of the cosipi, those that are > 1 
indicate angles that vanish on the segment of the caustic curve bordering 
the region, while those that are < — 1 indicate angles that become n on that 
segment. The pattern of 0's and 7r's uniquely identifies the region, as shown 
by Table GQ 

With these definitions, the imaginary part of the analytic continuation of 
$ PR is 

$pr = E J ^, (20) 

i 

where the sum runs over all six i. The quantity $pr vanishes on the caustic 
curve and becomes real and negative as we move into classically forbidden 
regions A or D, or real and positive as we move into regions B or C. The 6j 
symbol decays exponentially as we move into any classically forbidden region, 
a behavior that is captured by exp(— |$pr|) in all regions. Since $pr has only 
one sign in any of the four classically forbidden regions, its absolute value is a 
smooth function in those regions. Finally, the Ponzano-Regge approximation 



11 



in the classically forbidden regions is 



I 3l 32 312 

\ 33 H 323 

where the amplitude is given by 



A PR exp(-|$ PR |), (21) 



2Jl2n\V\ 



and where uqj is given by f|T8|) . The definitions we have made allow a sin- 
gle formula to be written down for all four regions, but it is easy to write 
four formulas for the four regions without the use of sign or absolute value 
functions. 



2.4 The Phase Space of the 6j -symbol 

The phase space of the 6j-symbol is the set of all closed figures that can 
be obtained by placing Jj, i = 1,2,3,4 end-to-end, for fixed values of the 
lengths of these vectors, modulo proper rotations. That is, the vectors must 
satisfy ([2]). This is the point of view of Kapovich and Millson[T6]. that 
recently has been further developed by Charles[13]. This space can also be 
derived by symplectic reduction [T7] from a model of four independent angular 
momenta built around Schwinger's[18] oscillators, much as in the treatment 
of Aquilanti et al[H] of the 3j-symbol. In our applications we will think 
of the lengths Ji, i = 1,2,3,4 as being given by the quantization condition 
([8]) in terms of the fixed quantum numbers ji, i = 1,2,3,4 appearing in a 
6j-symbol. For a given closed chain formed by J;, i = 1, 2, 3, 4, we can draw 
vectors Ji 2 and J 2 3 defined by ([3]) to obtain a tetrahedron. Thus the phase 
space can also be thought of as the set of all real tetrahedra, in which the 
lengths Ji, i = 1,2,3,4 are fixed. The lengths Jyi and J23, however, vary 
continuously between the limits (fTUl) . 

All such tetrahedra are generated if we let Ji 2 vary from Ji2, m in to Ji2, ma x, 
while for each value of J12 we let the dihedral angle 0i2 vary from — 7r to 
+7T. Here 0i2 is the interior dihedral angle, illustrated in Fig. that is 
uniquely defined in the interval — ix < <p\2 < ^ by requiring that <pi 2 = 0, n 
correspond to flat tetrahedra, and that < <pi2 < n correspond to tetrahedra 
with positive volume (this is the case illustrated in Fig. [5]). It is related 
to the exterior dihedral angle ipi2 used in the Ponzano-Regge formula by 
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1 012 1 + ^12 = 7T- The angle i2 distinguishes tetrahedra related by spatial 
inversion (i.e., time reversal), while ipi2 does not. Similarly, we could generate 
all real tetrahedra by varying J23 and the dihedral angle 023- The choice of 
J12, 012 for this process is arbitrary, but it gives us a definite convention for 
coordinates on the phase space of the 6j-symbol, namely, (Ji2,0i2)- 

The manifold of such tetrahedra, modulo proper rotations, is a sphere. 
To see this, define 

^12,avg = —(Jl2,ma,x + ^12,mm) 5 (23) 

and write 

K z = J\2 — t/l2,avg; (24) 

so that K z varies between —D/2 and +D/2 as J i2 goes from Ji2, m in to Ji2,max 
(see (TT2"j) ). Then define a polar angle 6 ) 12 by 

K z = (D/2)cos6 12 , (25) 

and set 

K x = (D/2) sin ^12 cos0i 2 , 

K y = (D/2) sin sin 0i2, (26) 

so that (K x , K y , K z ) are Cartesian coordinates on a sphere of radius D/2 
with spherical angles (6 I 12 , 0i2)- The azimuthal angle 0i2 on the sphere is the 
same as the interior dihedral angle in the tetrahedron. 

This is the 6j-sphere, on which the north pole is K z = D/2 or J12 = 
-7i2,ma X , the south pole is K z = —D/2 or Ji 2 = Ji2,min, and curves of constant 
J12 in general are small circles K z = const. It is illustrated in Fig. [61 which 
shows several curves of constant J 12 . Flat configurations correspond to 0i2 = 
or 7r, that is, they lie on the K x -K z plane (the great circle K y = 0). The 
hemisphere K y > (K y < 0) consists of tetrahedra of positive (negative) 
volume. Spatial inversion (i.e., time reversal) is a reflection in the plane 
K y = (it amounts to K y — > —K y ). 

Any quantity defined in terms of the tetrahedron that is invariant under 
proper rotations corresponds to a function on the 6j-sphere. For example, 
J23 is such a function, as is 023- Curves of constant J23 are illustrated in 
Fig. [71 The extrema of J23 are both flat configurations lying on the great 
circle K y = 0, with J23 = J23, m in on the semicircle K x > and J23 = J23,max 
on the semicircle K x < 0. This is apparent from figures such as diagram 
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X in Fig. HI which illustrates the case J23 = J23,min and which shows that 
012 — (not tt) at such a configuration. Two views of the 6j-sphere are 
given in Fig. [7] to show both extrema of J23, as well as curves of constant J23 
for intermediate values. Notice that the curves of constant J23 are not small 
circles. That is because the coordinates we are using on the 6j-sphere make 
the curves of constant J12 look simple (they are small circles), but not the 
curves of constant J 2 3- Had we based our coordinates on J 2 3 and 2 3 instead, 
the roles would be reversed. 

The diagrams in Fig. [7J were created in the following way. We set up a 
grid of coordinates ( J12, 4>u) by letting J i2 vary between Ji2, m in and Ji2,max, 
and for each value of J12, letting 12 vary between ±7r. For each value of 
J12 at 0i2 = 0, we set up the corresponding flat tetrahedron, such as the 
diagrams labelled Y — > X and X —>■ Z in Fig. HI both of which have 0i2 = 0. 
For other values of 0i2, we rotate the triangle 1-2-12 about the axis defined 
by J12 by angle 0i2, using the right-hand rule, while holding triangle 3-4-12 
fixed. This is the "butterfly" motion of the tetrahedron associated with the 
axis Ji 2 . This creates a tetrahedron of the desired dihedral angle 0i 2 . Then 
we compute J23 for that tetrahedron. In this way, we set up an array of J 2 3 
values on the grid. Finally, we draw the contour lines of J23 for this grid, 
and plot them on the surface of the sphere. 

There are two ways to compute Poisson brackets on the 6j-sphere. First, 
let F and G be two functions of the J,, i = 1,2,3,4. Then the Poisson 
bracket is the usual one in classical mechanics for a set of independent angular 
momenta, 




For example, any function of Ji and J2 has vanishing Poisson bracket with 
any function of J3 and J4. 

Hamilton's equations can be expressed in terms of Poisson brackets. Let 
if be a Hamiltonian with evolution parameter (the "time") A, and let F be 
any function of J,, i = 1,2,3,4. Then the rate of change of F along the 
orbits of H is 

% = {F,H } . (28) 
For example, if we take H = J 12 = |Ji + J2I and F = some component of 
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one of the Jj, we find 



dJi _ j eu x Ji, i = 1,2, , s 

dX \ 0, i = 3,4, l2Jj 

where ei2 is the unit vector in the direction J 12. The A-evolution is a rotation 
of vectors Ji and J2 about the axis e±2, while J3 and J 4 remain fixed. This 
is the "butterfly" motion mentioned above, and A is the angle. If the initial 
conditions are chosen so that A = when i2 = 0, then A = 0i2- On the 
6j-sphere the orbits of Jyi are the curves Jyi = const, the small circles seen 
in Fig. EJ 

Similarly, J23 generates another butterfly motion, in which vectors J2 
and J3 rotate about the axis e23 by the right-hand rule, with dihedral angle 
023 as the parameter of evolution. Vectors Ji and J4 remain fixed during 
this motion. The orbits on the 6j-sphere are curves of constant J23, some 
examples of which are illustrated in Fig. UJ 

Since a Hamiltonian and its evolution parameter are canonically conju- 
gate variables, it follows that ( J12, $12) are canonically conjugate variables 
on the 6j-sphere. Since J 12 differs from K z by a constant, we can equally 
well use (K z , 0i 2 ) . Thus another way to compute the Poisson bracket of any 
two functions defined on the 6j-sphere is 

{F,G} = ^^-^^. (30) 

(9012 dJ U 8Ji2 <90i2 ' 

This only applies to rotationally invariant functions, while ( l27j) can be used 
for any functions of the Jj, i = l,2, 3, 4. Also, expressing functions in terms 
of J12 and 0i2 is usually difficult, so in practice (127)) is more useful than (13"0)) . 

But (130)) does show that the area of a closed curve on the 6j-sphere can 
be computed as 

Area = f J 12 d(f) 12 , (31) 



with due attention to the singularities of the (Ji2,0i2) coordinates (there 
is no global symplectic 1-form on the sphere). The area is D/2 times the 
solid angle subtended by the closed loop. This is not the area on the surface 
of a sphere of radius D/2 in (K x , K y , K z ) space, computed by Euclidean 
geometry, which would be {D/2) 2 times the solid angle, but it is the correct 
measure of area from the standpoint of semiclassical mechanics. 

In particular, the total area of the sphere is (Air) (D/2) = 2ttD, or D 
Planck cells of area 2tt each (it would be 2irh in ordinary units). This is 
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what we expect for a semiclassical phase space representing the subspace Z 
of the Hilbert space of four angular momenta, which contains D quantum 
states. 

Moreover, the Bohr-Sommerfeld rules say that the quantized values of 
J12 should be given by orbits Ji 2 = const whose area is (n + |)(27r). On the 
sphere there is no way to distinguish the interior and the exterior of a loop, 
but the Bohr-Sommerfeld rule is the same in either case since the total area of 
the sphere is an integer times 2tt. Since the small circle Ji 2 = const encloses 
area (27r)(Ji 2 — </i 2 , m in) about the south pole, the quantized orbits are those 
for which J i2 = Ji2,min + \ + integer. By ([8]) this gives the exact quantized 
values of ji 2 , as indicated by (jlj). The minimum and maximum quantized 
values of Ji 2 are one half unit away from the values at the south and north 
poles, respectively, corresponding to the one-half unit margin between the 
quantized spots in Fig. [3] and the bounding values of Ji 2 . Figure [6] illustrates 
the quantized orbits of J i2 for the case D — 5, numbered to 4 as Ji 2 
increases (orbit lies close to the south pole and cannot be seen in the 
figure). 

Similarly, the quantized orbits of J23 are those satisfying J 23 = J 2 3 im i n + 
\ + integer. These are labeled to 4 in Fig. [7J for the same case [D = 
5) illustrated in Fig. [6j The dihedral angle 2 3 is an angle parameterizing 
position along the curves J 23 = const, although it is not an azimuthal angle 
in (K x , K y , K z ) space. 

But the logic that leads to the conclusion that ( Ji 2 , 0i 2 ) are canonical 
coordinates on the sphere applies also to (J 2 3, </> 2 3), so there is a canonical 
transformation connecting (Ji 2 ,0i 2 ) and (J 2 3,0 2 3). According to Miller's 
theory [19], the i^-type generating function of this canonical transforma- 
tion is the phase of the semiclassical matrix element O12 1 J23) - (We follow 
Goldstein's [20J conventions for classifying generating functions.) 

Miller's theory leads to difficult integrals in cases like this (it certainly 
does for the 3j-symbol), and it has never been carried through for the 6j- 
symbol, as far as we know. But it is certain that the i^-type generating 
function that would result would be the Ponzano-Regge phase $pr, to within 
an additive constant. Moreover, Miller's theory shows that the amplitude 
determinant in the Ponzano-Regge formula is given by 

1/2 

(32) 



.4 



PR 



PR 



dJ 12 dJ 23 

to within a multiplicative constant. In this formula, the lengths Jj, i 
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1, 2, 3, 4 are considered fixed, and the dihedral angles ipi that appear in (fTST) 
are considered functions of all six lengths Jj. 

The Ponzano-Regge amplitude was first derived by Wiener [21] . who had 
the intuition that the probability in making a measurement of j'23 for given 
value of j'12 should be uniformly distributed in the angle (pu. That this is so 
follows from standard semiclassical theory and the fact that 4> 12 is conjugate 
to J\2- This amplitude is also inversely proportional to the square root of a 
Poisson bracket, 

fe,^} = J '-' J2 , Xj3 ' = 7 ^ r . (33) 

J\2 ^23 ^12^23 

which was computed using (12T1) . The volume factor V is the part of the am- 
plitude that was obtained by Wigner, while the factor J12J23, when replaced 
by (ji2 + |)(j23 + |), is needed to convert from the unitary matrix element 
(712 1.723) to the 6j-symbol (see Edmonds pL J Eq. (6.2.10)). The use of Poisson 
brackets for computing amplitude determinants is discussed in Littlejohn[22] 
and in Aquilanti et al[14j. 

The relation between Fig. [3] and the phase space of the 6j-symbol is the 
following. A point inside the square of Fig. [3] specifies values of J± 2 and J 2 3 
that are allowed by the inequalities (Q. These in turn specify two curves 
in the phase space, one of constant J\ 2 and the other of constant J 23 . If 
these curves intersect, as in part (a) of Fig. El then we are in the classically 
allowed region. In that case, the intersection points, labeled P and Q in 
the figure, are the stationary phase points of the semiclassical evaluation 
of the matrix element (J12 1 J23) - These points represent two tetrahedra that 
are mirror images of each other (they are related by time reversal). The 
total semiclassical matrix element (the Ponzano-Regge formula) is a sum of 
contributions from these two tetrahedra, which are complex conjugates of 
each other. Thus, the semiclassical matrix element is real. 

If the two curves do not intersect, as in part (b) of Fig. [8], then we are in the 
classically forbidden region. Both curves are manifolds of real tetrahedra, one 
with a fixed value of J± 2 , the other of J23, but since they do not intersect, there 
is no real tetrahedron that simultaneously has both given values of J\ 2 and 
J23. In this case the analytic continuations of the curves into complex phase 
space (a complexified sphere) do intersect, and these intersections represent 
the stationary phase points in the classically forbidden region. We make no 
attempt to sketch the complexified phase space, however. 
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In addition to its interpretation as a generating function, the phase of 
semiclassical matrix elements such as (ji2 IJ23) is geometrically one half the 
area enclosed by the intersection of the quantized classical orbits in phase 
space[22]. (The relative phase between the two branches of the WKB solution 
is the area, but this is shared between two exponentials to create a cosine 
term. Thus the argument of the cosine is one half the area.) In part (a) of 
Fig. [S] the area in question is the shaded region (a "lune"). 

If the two curves J12 = const and J23 = const are tangent, then we are 
at a caustic. A caustic implies a flat tetrahedron of zero volume, so such 
tangencies can occur only in the plane K y = 0. 

The different types of caustics than can occur are illustrated in Fig. O 
Part A of Fig. is obtained from the curves of Fig. [H] by adjusting the 
J12 or J 2 3 values to create a tangency. As expected, it lies between the 
classically allowed region (part (a) of Fig. [HJ and the classically forbidden 
region (part (b) of Fig. [8]). The point of tangency T is the caustic point. 
From part A of Fig. [9] we move into the classically allowed region if we either 
decrease J i2 or increase J23. Thus we see that it corresponds to region A of 
Fig. El 

If we allow J12 in part A of Fig. [9] to decrease, the small circles Ju = 
const sweep down from the north pole through the oval J23 = const, passing 
through the classically allowed region, until a tangency is reached at the 
lower point of the oval J23 = const, where the small circle J 12 = const is 
close to the equator. This is another caustic, illustrated in part B of Fig. 
The shaded area is the continuation of the shaded area in part (a) of Fig. [HJ 
and again T is the caustic point. From this caustic, we pass back into the 
classically allowed region if either J12 increases or t/23 increases, so we are in 
region B of Fig. [3J 

If from part A of Fig. M we allow J23 to increase, then the curves J23 = 
const sweep through the small circle Jyi = const about the north pole, finally 
reaching another tangency on the other side where K x < 0. The result is 
illustrated in part C of Fig. [9], where again point T is the caustic point. From 
this caustic we pass back into the classically allowed region if we let either 
J12 or J23 decrease, so this corresponds to region C of Fig. [31 

As J23 increases from its value in part C of Fig. [91 the curve J23 = const 
shrinks down around the point J 2 3 = ^23,max on the semicircle K y = 0, 
K x < 0. Then allowing J12 to decrease, the small circle around the north 
pole moves south, passing through the curve J 2 3 = const, producing finally 
a tangency T on the other side, as illustrated in part D of Fig. Now the 
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shaded area (twice the Ponzano-Regge phase plus a constant) covers nearly 
the entire sphere. From this configuration we pass back into a classically 
allowed region if either J23 decreases or J12 increases, so we are in region D 
in Fig. [31 

In the case of an ordinary oscillator with a flat phase space (the plane), 
the difference in the actions between the two turning points is one half the 
area of the orbit, and has the form (n + |)7T, where n is an integer. As 
explained in the introduction, this is a requirement for the existence of a 
uniform approximation of the Weber function type. The analogous statement 
for the 6j-symbol with the spherical phase space is sometimes true, and 
sometimes not. A case where it is true is obtained from diagrams A and 
B of Fig. in which we regard j'23 fixed and j'12 variable. As J\i decreases 
from the north pole (its maximum value), we first encounter a caustic of 
the type A, where the area of the lune is zero. Continuing to decrease Ji 2 , 
we pass through the classically allowed region, finally encountering a caustic 
of the type B (the lower turning point), where the area of the lune is the 
quantized area of the oval J23 = const. This area has the form (n + |)27r, 
so the differences in the actions at the two turning points is quantized. This 
implies that the difference in the Ponzano-Regge phases $pr between the 
two turning points is also quantized. 

A case where the differences in the actions is not quantized and a uniform 
approximation of the Weber function type does not exist is obtained when 
J23 has a value such as that illustrated in part C of Fig. O In this case, 
as we let J12 decrease from its maximum value at the north pole the first 
caustic we encounter is of type C, where the area of the lune is the area of 
the curve J12 = const (the shaded area in part C of the figure). This area 
is not quantized, since the value of J i2 at a caustic is not quantized. As J\2 
decreases, we eventually reach the lower caustic of type B, where the area of 
the lune is the quantized area of the orbit j 2 3 = const. Thus, the differences 
between the areas, one quantized, the other not, is not quantized. A case 
like this (with caustics of the type B and C) was illustrated in Fig. HJ 

We will now show that the <i-matrices have a phase space and an orbit 
and caustic structure that are identical, from a topological standpoint, to 
those of the 6j-symbol. 
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3 The d-matrices 



3.1 Quantum Mechanics of the d-Matrices 

The d-matrices are defined by 

di m ,((3)^{m\U y ((3)\m'), (34) 

where U y (/3) = exp(—ij3J y ) is a rotation operator with Euler angle (3 about 
the y-axis, and \m) and \m') are standard rotation basis states (eigenstates 
of J z ). To indicate both the operator and the quantum number, we will write 
these states as | J z : m) and \J Z : m'). By conjugation the rotation operator 
U y (/3) rotates the angular momentum vector, 

Uy{(3)^Uy{(3) = Ry(P)J, (35) 

where R y {f3) is the 3x3 rotation matrix for an active rotation about the 
y-axis. We define 

/ cos (3 sin (3 \ f \ / cos (3 \ 
n = R y ((3)z =\ 1 = , (36) 

\ - sin/3 cos/3 / \ 1 / \ sin/? / 

as illustrated in Fig. [TO], so that 

(n ■ 3)U y ((3)\J z : rri) = U y ((3)(z ■ J)| J z : m') = m'U y ((3)\J z : m'), (37) 

where we use ( 135]) and 

n • [R y ((3)3] = [Ryi^ii} ■ J = z • J = J z . (38) 

Therefore U y (f3)\J z : m') is an eigenstate of n • J = J n with eigenvalue m', 
and we will write 

U y (f3)\J z :m') = \J n :m'), (39) 

so that 

d? mm ,{(3) = {J z :m\J n :m!). (40) 

In this way the (i-matrix is written as a unitary matrix element connecting 
the eigenstates of two different operators. This is the starting point for 
Miller's {19] theory of semiclassical matrix elements, as well as our own|221[T4"] 
treatments of the same subject. 
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3.2 Classical and Semiclassical Mechanics of the d-ma- 
trices 



References on the semiclassical approximation for the d-matrices include 
Brussaard and Tolhoek[23j, Ponzano and ReggeJTU], Braun et al[21] and So- 
kolovski and Connor [25]. In the following we emphasize geometrical aspects 
of the problem not covered by these authors. 

The classical phase space for d 3 mm ,{(5) is a sphere ("the d-sphere") in 
angular momentum space of radius |J| = J, where 

J = 3 + \- ( 41 ) 
The area of a loop on the surface of the sphere is given by 

Area = f J z d<f>, (42) 



where is the azimuthal angle, again with due consideration of the singu- 
larities of the coordinates (<p,J z ). That is, if the loop subtends solid angle 
Q, then the area is Jfl (not J 2 Q, as in Euclidean geometry). The total area 
of the sphere is therefore An J = (2j + l)(27r), that is, the sphere consists 
of 2j + 1 Planck cells, corresponding to the 2j + 1 basis states \ J Z : m) or 
\J n : m!). Curves of constant z • J = J z and n • J = J n are small circles 
centered on the axes z and n, respectively, as illustrated in Fig. [TT1 

The Poisson bracket of two functions F and G of J is (J27j) with a single 
term in the sum, 



dJ dJ 



or, equivalently, for functions of (0, J z 



For example, using (143]) we find that Hamilton's equations for Hamiltonian 
J z = z • J with evolution parameter A are 

^ XJ . (45) 

The motion is a rotation about the z-axis, so the orbits are the small circles 
J z = const, as expected. The parameter of the orbit is A = 0, so and J z are 
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conjugate variables, as indicated in (j44p . Similarly, J n generates rotations 
about the axis n. 

Let (9, 0) be the usual spherical coordinates referred to the axis z, and 
let (9', 0') be an alternative set referred to the axis n. That is, the (#', 0') 
coordinates of a point (x, y, z) are the same as the (9, 0) coordinates of the 
inverse rotated point R y (/3)~ 1 (x,y, z). Thus the coordinate transformation 
(9, 0) -> (0', 0') is specified by 

sin 9 cos = cos /3 sin cos 0' + sin /3 cos 0', 
sin 9 sin = sin 9' sin 0', 

cos 9 = — sin /3 sin 6 1 ' cos 0' + cos /5 cos 0'. (46) 

The azimuthal angle 0' is conjugate to J n , so both (0, J z ) and (0', J n ) are 
canonical coordinates on the sphere. The i^-type generating function of 
the canonical transformation between these coordinates is the phase of the 
semiclassical approximation to the (i-matrices, according to Miller's theory. 
This aspect of the problem has been developed by Sokolovski and Connor |25j. 

The classical observables J z and J n are functions on the <i-sphere that 
vary continuously between the limits, 

~J<Jz,Jn< +J- (47) 

The quantized orbits of J z and J n are those enclosing n+| Planck cells where 
n is an integer. This implies J z = m and J n = ml with the usual rules for 
quantum numbers m and m', 

— j < m,m' < +j, (48) 

in integer steps. Thus the maximum and minimum values of m and ml lie 
one half unit away from the maximum and minimum values of the classical 
observables J z and J n , as illustrated in Fig. [12j This figure may be compared 
to Fig. [3] for the 6j-symbol. See also Fig. 1 of Braun et al|24j. 

When the J z -orbit and the J n -orbit intersect one another as in part (a) 
of Fig. [H], then we are in the classically allowed region of the ^-matrices. 
There are generically two intersection points related by a reflection in the 
plane J y = 0, marked by unit vectors from the origin a and a' in Fig. [T3j 
We concentrate on intersection a, for which J y > 0; at the other intersection 
a' we have J y < 0. In the coordinate systems (6*,0), (6 )/ ,0 / ), the 9 and 9' 
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coordinates of intersection a are given by 



cosO = — = r , cos6> = — = r , (49) 

J j+\ J 3+\ 

where J z and J n label the two small circles and where the second form applies 
if J, J z and J n take on their quantized values. As for the and <$' coordinates 
of intersection point a, they can be obtained by solving PET) , assuming 9, 9' 
and f3 are given. This gives 

cos 9' — cos j3 cos 9 cos 9 — cos /3 cos 9' 

cos0= . . , cosr/ = . , , 50 

sin p sm # sin p sin 

where we write r] = n — 4>' as illustrated in Fig. [131 Equations fl50|) uniquely 
determine <fi and ^ in the interval [0, 7r] (assuming that the J z - and J n -orbits 
actually intersect). 

Figure d3] draws attention to the spherical triangle defined by z, n and 
a, whose sides are arcs of great circles subtending angles 9, 9' and (3. Equa- 
tions (|50|) are the law of cosines for spherical triangles, applied to the interior 
angles and 77, as shown in the figure. As for the third interior angle, we 
define k as the opening angle of the lune (the shaded area), as illustrated 
in the figure. Then it is easy to show that the third interior angle of the 
spherical triangle at vertex a is 7r — k. For this angle the law of cosines gives 

. . cos (5 — cos 9 cos 9' . 

COS(7T — K) = = — COS K, (51) 

sm 9 sm 9' 

determining k also uniquely in the interval [0, 7r]. 

We define $^ as one half of the area of the shaded lune seen in Fig. [TSl 
which is also the i^-type generating function of the transformation (9, 0) — > 
(9',(j)') (Sokolovski and Connor [25J). Then we have 

$d = J(k — <j)cos9 — i]cos9') 

= Jk - J z 4> - J n r) = (j + -)k- m<f) -m'r), (52) 

where the final form applies if J, J z and J n take on their quantized val- 
ues. The result is the sum of angular momentum quantum numbers times 
dihedral angles, that is, the interior angles of the spherical triangle formed 
by (z, n, a) are also the interior dihedral angles of the tetrahedron or par- 
allelepiped formed by those vectors. The analogy with the Ponzano-Regge 
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formula is more transparent if we use the exterior dihedral angles k, it — <p 
and 7r — 77, so that 

1 

= (J + + m ( 7r ~ 0) + 171 fa ~ V) ~ ( m + rn')'^, (53) 

in which the first three terms look like the sum (TL5I) . while the final term just 
produces a phase factor (— l) m+m in the asymptotic formula. 

It is straightforward to prove (152j) by elementary geometry, but another 
proof, based on symplectic reduction of Schwinger's oscillator model of an- 
gular momentum[18] (essentially the Hopf fibration), leads to the following 
generalization. Let a polygon on the unit sphere be specified by vertices 
(vi, . . . , v n ) connected by arcs of small circles, where the small circle pro- 
ceeding from Vj to v i+1 is obtained by rotating Vj about axis n, by angle (pi, 
using the right-hand rule. Also, let Ki be the interior angle between the two 
small circles meeting at Vj. Then the solid angle of the interior of the poly- 
hedron, defined as the region to the left as we move along the small circles, 
is 

ft = 27r-f;[(7r-K < ) + (v i .nO&]. (54) 

i=l 

Special cases of this formula include the solid angle of a spherical triangle 
(with sides that are great circles), Q = K\ + + K3 — it, and (|52|) . for which 
n = 2, K\ = K2 = k, and $^ = Q/2. Equation (j54j) can also be derived as a 
special case of the Gauss-Bonnet theorem. 

The spherical triangle formed by (z, n, a) plays another role. We define 
Vd as the volume of the parallelepiped spanned by these three vectors, which 
can be written in a variety of ways, 

Vd = (z x n) ■ a = sin (3 sin 6 sin <p 

= sin j3 sin 0' sin tj = sin 6 sin 8' sin /t, (55) 

where we use the law of sines for the final three expressions. One of these 
equalities is equivalent to the ^-component of (H6]l . The square of Vd is the 
determinant of the Gram matrix formed by vectors (z, n, a), 

(z • z z • n z • a \ / 1 cos (3 cos 9 

n • z n ■ n n ■ a = det cos (3 1 cos 9' 
a-z a-fi a-ay y cos 9 cos 9' 1 

= l + 2cos/3cosflcosfl'-cos 2 /3-cos 2 fl-cos 2 fl'. (56) 
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Since we are working with the volume of the parallelepiped instead of the 
volume of the tetrahedron, there is no factor of 6 in (j55p or of 6 2 = 36 in 
(156]) . The volume Vd appears in the amplitude of the asymptotic formulas 
(J57D and (jMj). 

Including all the details (Maslov indices, phase conventions, etc), the 
asymptotic expression for the ci-matrix in the classically allowed region is 

Cn'(/?) = A*cos($ d -£), (57) 

where the amplitude is 

A d = , l ] (58) 

According to Ref. [22] the amplitude of the WKB approximation for the 
d-matrix should be proportional to the inverse square root of the Poisson 
bracket { J z , J n }, evaluated at the intersection between the two orbits, J = 
J a. Indeed, using fj43|) . we see that it is: 

{J z ,J n } = (z x n) • J = JV d = sin pj y 

= J sin (3 sin 6 sin <p. (59) 

The caustics of the <i-matrices occur when the small circles J z = const 
and J n = const are tangent, or, equivalently, when the vectors (z, n, a) are 
linearly dependent so that = 0. Multiplying (1561) by J 2 , using (H9|) and 
setting the result to zero gives us the equation of the caustic in J z -J n space, 

J 2 Z + Jl - 2 J z J n cos (3 - J 2 sin 2 (5 = 0, (60) 

an ellipse whose axes are oriented 45° to the J z -J n axes, and whose semimajor 
and semiminor axes are a/2 cos(/9/2), \/2 sin(/3/2). An example is illustrated 
in Fig. [TJl see also Fig. 1 of Ref. [21]. The ellipse touches the boundary 
defined by the classical limits (T4T1) at four points, creating four classically 
forbidden regions labeled ABCD in Fig [T2J The square of the volume V d 2 is 
negative in the classically forbidden regions, and itself is imaginary there. 

Another point of view on the caustics is to hold J z and J n fixed, thereby 
fixing the sizes of the two small circles, and to vary (3, which moves the 
position of the small circle J n = const. Then the small circles are tangent at 
the turning points f3 = (3\ or fa, where < (3\ < (32 < vr, and where 

j3 l = \e-9'\ 1 f3 2 = mm(6 + 6',2n-6-6'). (61) 
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The classically allowed region is (3i < (3 < (3 2 , while the two classically 
forbidden regions are < /3 < (3% and (3 2 < (3 < tt. 

The four types of tangencies of the two small circles are illustrated in 
Fig. [T4] In all four parts of the figure, T is the caustic point (the point of 
tangency). In part A we are at the upper turning point (3 = [3 2 , because if (3 
decreases we obtain two intersection points and are in the classically allowed 
region. In fact, this is the case j3 2 = 9 + 9' < tt. Or if we hold (3 fixed 
but decrease either J z or J n , again we enter the classically allowed region, 
since one or the other of the two small circles expands and the tangency 
develops into two intersection points. Thus part A of Fig. [TH corresponds to 
the corner A of Fig. [T2J In part B of Fig. [TH we are at the lower turning 
point P = Pi = 9 — 9' > 0, since if (3 increases we move into the classically 
allowed region. The same happens if we hold (3 fixed and either increase J z 
or decrease J n , so this corresponds to corner B of Fig. [12] In part C of Fig. [HI 
we are at the lower turning point j3 = (3\ = 9' — 9 > 0, which corresponds to 
corner C of Fig. [12] since we enter the classically allowed region if either J n 
increases or J z decreases. Finally, in part D of Fig. [TJ] we are at the upper 
turning point (3 2 = 2tt — 9 — 9' < tt, which corresponds to corner D of Fig. [T2] 
since we enter the classically allowed region if either J z or J n increases. 

The four types of tangencies of orbits for the 6j-symbol, illustrated in 
Fig. [9] are topologically identical to the four types for the <i-matrices, il- 
lustrated in Fig. dH Similarly, the four classically forbidden regions of the 
6j-symbol, illustrated in Fig. [3] are in one-to-one correspondence with the 
four classically forbidden regions of the <i-matrices, illustrated in Fig. [12] 
Comparing Figs. [3] and [121, we see that the labelings of the corners by ABCD 
are not the same; but this is because the point on the <i-sphere of maximum 
J n , namely, the point in the direction n, corresponds to the point on the 
6j-sphere of minimum J 23 . If the J23 axis in Fig. [3] had been drawn increas- 
ing downward instead of upward, the labels on all four corners (classically 
forbidden regions) of both Fig. [3] and Fig. [12] would coincide. 

Referring to Fig. [12], if we hold (3 fixed and vary J z or J n , moving from 
the interior of the ellipse (the classically allowed region) to the boundary (the 
caustic), then all the angles 0, r] and k approach either or tt, depending 
on which segment ABCD of the boundary (the caustic curve) we approach. 
The values of these angles on the caustics are summarized in Table [2] For 
uniformity of notation, we write at, i = 1,2,3 for k, cf) and rj, as indicated 
in the table, and similarly we write ki, z = 1,2,3 for j, —m, —m', where the 
signs are the same as in the three terms of the expression flo^]) for Also 
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Table 2: Values of the angles k, <ft and rj on caustics of type ABCD, also 
integer i/j, for four caustic types. 

shown in the table is the integer v& for the four classically forbidden regions, 
defined by 

Vd = ^h = integer, (62) 

i 

where the sum is only taken over i such that oti = it. The definition is similar 
to that of in (ITS1) . and used for a similar purpose, that is, the asymptotic 
form of the (^-matrices in the classically forbidden regions carries a phase 
(— l) Ud , effectively due to the analytic continuation of 

The angles oti are extended into the classically forbidden region in a man- 
ner exactly like that used for the ipi in the case of the 6j-symbol, as explained 
below (fl8l) . That is, if a« = on the segment of the caustic curve adja- 
cent to a given classically forbidden region, then we define aj = ia,i, where 
Q-i = cosh - (cosctj) is real and positive; while if ol{ = tt on the caustic curve, 
then we define an = tt + ion, so that oci = — cosh -1 (— cos ctj) is real and 
negative. In the classically forbidden regions, the quantities cosaj, given by 
f l50|) and fl5T|) . lie outside the interval [—1, —1]. 

We now define a quantity related to the analytic continuation of $^ into 
the classically forbidden regions, 

3 3 

= ^2 kiOCi = ^2ki sign(cos Q!j) cosh -1 | cos o;j| . (63) 

i=l i=l 

In spite of the absolute value and sign functions, 5^ is smooth over any given 
classically forbidden region. This is important for root finders that rely on 
smoothness, such as the Newton-Raphson method. The quantity is zero 
on the caustic boundary, and real and negative as we move into classically 
forbidden regions B and C, and real and positive as we move into classically 
forbidden regions A and D. These are the same rules as for $pr. The d- 
matrices decrease exponentially as we move into any classically forbidden 
region, a behavior that is captured by exp(— \<&d\) i n & h cases. 
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Finally, the asymptotic expression for the <i-matrices in the classically 
forbidden region is 

d Lm> = A d exp(-\$ d \), 

where 

2^J(7r/2)J\V d \' 

4 The Uniform Approximation 

4.1 Remarks on Uniform Semiclassical Approximations 

The traditional method of constructing uniform semiclassical approxima- 
tions, the "method of comparison equations," is reviewed by Berry and 
Mount [26], with citations to earlier literature. In this method one takes 
a one-dimensional Schrodinger equation (a second-order differential equa- 
tion in x) and performs a coordinate transformation X = X(x) to create 
a new Schrodinger equation in X which, after the neglect of terms of order 
h 2 , becomes a standard, solvable equation. The most common standard or 
"comparison" equations in practice are the differential equations for Airy or 
Weber (parabolic cylinder) functions. 

Since both the 6j-symbols and the <i-matrices satisfy second-order differ- 
ence equations, it is likely that a kind of discrete version of the method of 
comparison equations, along the lines of the discrete WKB theory used by 
Schulten and Gordon [TT] and refined by Braun[15j, could be used to construct 
a uniform approximation for the 6j-symbol. We have not constructed our 
approximation in this way, however, and if we had, it is likely that we would 
have missed much of the geometry discussed above. Also, that approach 
would have produced a uniform approximation for the 6j-symbol only for a 
fixed value of ju or j'23, not over the whole range of both variables as we have 
done here (it would not have produced a uniform approximation in terms of 
o?-matrices) . 

Hiding slightly beneath the surface of the usual method of comparison 
equations is a transformation between the phase spaces of the original prob- 
lem and the standard problem. This transformation is X = X(x) and 
P = (dx/dX)p, where the first part is the coordinate transformation used 
in the method (a "point transformation") and the second part is the usual 
lift of a point transformation into a canonical transformation. The second or 
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momentum equation can also be written pdx = P dX, which by integration 
gives 



where s and S are the actions of the original problem and the standard 
problem, respectively. In fact, this equation (the equality of the actions) 
specifies the coordinate transformation X = X(x). The geometry is illus- 
trated in Fig. [T5J in which the method of comparison equations is used to 
map a quantized curve of a nonlinear oscillator (the Morse oscillator, part 
(a) of the figure) into a quantized curve of a standard problem (the harmonic 
oscillator, part (b) of the figure). Equation ([66]) implies the equality of the 
shaded areas in the figure, which in turn determines the function X(x). In 
the figure, X = X(x ). 

The function X(x) can be analytically continued to the classically forbid- 
den regions where p, P, s and S all become complex, but the transformation 
X = X(x) is real and in fact forms a single, smooth (usually analytic) co- 
ordinate transformation throughout both classically allowed and forbidden 
regions. See Fig. 2 of Ref. [27J for a plot of the function X(x) in one example, 
which shows its completely smooth behavior as one passes from classically 
allowed to classically forbidden regions. The equation s(x) = S(X) is ob- 
vious in a sense: actions are areas, and areas are preserved by canonical 
transformations. 

We have been interested in the generalization of the method of compari- 
son equations to a class of canonical transformations that is larger than the 
point transformations. Since canonical transformations are the semiclassical 
representatives of unitary transformations, the idea is to carry out a unitary 
transformation on the original problem such that the transformed problem 
has a standard form, to within errors of order h 2 , where the choice of unitary 
transformation is guided by geometrical criteria in the classical phase space. 
The standard problem is sometimes referred to as a "quantum normal form." 
Examples of quantum normal form calculations and an interesting perspec- 
tive on Bohr- Sommerf eld or torus quantization may be found in Refs. [28, 29J. 
In those references only the problem of determining eigenvalues is considered, 
but in the present application we are interested in the transformation of the 
wave functions, a problem that involves extra features. As it turns out, it 
is not only necessary to carry out a unitary transformation (which at the 
classical level maps a pair of orbits into another pair that are in standard 
form), but also a certain nonunitary transformation (to make the densities 




(66) 



29 



of probability on the orbits come out in standard form). 

In the usual method of comparison equations, the uniform approximation 
for the exact solution ip(x) is given by 

Ms) « ( 67 ) 

where ^(X) is the standard solution of the standard problem, and a(x) and 
A(X) are the amplitudes of the two semiclassical approximations, a(x) for 
the original problem and A(X) for the standard problem. Both amplitudes 
diverge at the caustics, but their ratio has a definite limit and in fact is 
smooth everywhere across both classically allowed and forbidden regions. 

Similarly, it turns out that the uniform approximation for the 6j'-symbol 
in terms of <i-matrices is given by 

If f f } - (-ir*^0(/3), (68) 

where A PR is given by ffTTl) or ff22l) and A d by fl58l) or fl65l) in the classically 
allowed or forbidden regions, respectively, and where (— l)" 6 * is an extra 
phase defined in (17^1) . In other words, the uniform approximation for the 6j- 
symbol is of the same form (16"T|) that emerges from the method of comparison 
equations, apart from an extra phase. 



4.2 The Details of the Uniform Approximation 

At the heart of the new uniform approximation is a smooth, area-preserving 
map between the 6j-sphere and the ci-sphere that is parameterized by fixed, 
quantized values of Jn and J23, related to quantum numbers j'12 and j'23 by 
([8]). The area of the 6j-sphere and that of the ci-sphere must be equal, which 
implies 

D = 2j + 1, (69) 

where D is given by (J7J) or (IT2|) . Thus the value of the parameter j of the 
d-matrix in (168]) is determined. 

The quantum numbers j'12 and j'23 determine a specific pair of orbits on 
the 6j-sphere, the small circle J12 = jn + §> and the oval J 2 3 = J23 + |- The 
map is required to map the small circle Ji 2 = J12 + |, a quantized orbit, onto 
a small circle J z = const on the rf-sphere. Because area is preserved, the 
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small circle on the rf-sphere must also be quantized, and contain the same 
area about the north pole as the original small circle on the 6j-sphere. That 
is, we must have j 12 , max - ji2 = j - m, or, 

m = ju ~ jl2,avg, (70) 

where ji2,av g = C/i2,min + ji2,max)/2. This determines the quantum number m 
in fl68l) . We only require our map to map this specific small circle of constant 
J12 on the 6j-sphere onto the corresponding small circle on the (i-sphere; 
other small circles of constant J 12 on the 6j-sphere, for other values of J12, 
are not mapped to small circles of constant J z on the (i-sphere. 

Similarly, we require the map to map the quantized oval J23 = J23 + \ on 
the 6j-sphere onto a small circle on the (i-sphere that is centered about some 
direction n that lies on the semicircle J y = 0, J x > on the (i-sphere. The 
direction n is a function of the angle /3, which will be specified momentarily. 
Because area is preserved, the new small circle on the (i-sphere will be a 
quantized orbit J n = m', enclosing the same area about the axis n as the 
oval J 2 3 = J23 + I encloses about the point J 23 = J23,min- Since the minimum 
of J23 corresponds to the maximum of J n , the quantum number m' satisfies 

323 - J23,min = j - m' ', OT, 

m 1 = j 2 3,avg - 323- (71) 

This determines the quantum number m' in (168|) . We only require our map to 
map this specific oval of constant J23 on the 6j-sphere onto the corresponding 
small circle on the ci-sphere; other ovals of constant J23 on the 6j-sphere, for 
other values of J23, are not mapped to small circles of constant J z on the 
(i-sphere. 

Finally, the parameter f3 is determined by requiring that the area of the 
lune on the 6j-sphere should equal the area of the lune on the (i-sphere. 
Effectively, we rotate the small circle J n = wl until the two areas are equal. 
This is in the classically allowed region; in the classically forbidden region, 
the analytic continuations of the areas on the two spheres are set equal. In 
this way, if the 6j-symbol is in the classically allowed region, then so is the 
(i-matrix, and vice versa. This condition is the analog of f|66|) in the standard 
method of comparison equations. 

We take the classically allowed region first. We shall show elsewhere that 
the Ponzano-Regge phase is related to the area of the lune on the 6j-sphere 
by 

$ PR = -(Area of lune) + $ , (72) 
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where $o is an extra phase that is related to the topology of loops in the 
phase space of Schwinger's oscillators, as discussed in Ref. [H]. Without 
going into this, we can determine <£>o by evaluating both the area of the lune 
and $pr at any point in the classically allowed region or on the caustic curve, 
as in Fig. [31 since $pr is a continuous function of position inside and on that 
curve. A point on segment A of the caustic boundary is convenient, since the 
area of the lune vanishes there, as shown by part A of Fig. O In this way we 
find 

3 

$0 = (Jl + -h + ^3 + Ja + Jl2 ~ </l2,max)7T = (>cx + ~)^, (73) 

where 

^ex = jl + 32 + h + k + 312 - jl2,max- (74) 

Note that z/ ex is an integer. As a check one can evaluate $pr and the area of 
the lune at other points on the caustic boundary (segments BCD), and see 
that the answer for $ agrees with (|73|) . 

As for the area of the lune on the ci-sphere, it is a function of /3 and is 
twice $d, given by (1521 . Altogether, the equation that must be solved for (3 
in the classically allowed region is 

$pr = $ d (/3) + $o- (75) 

In spite of the complications arising from the extra term <3>o> the geometrical 
meaning of (175|) is simple: the areas of the lunes on the two spheres are equal. 
Taking cosines, (1751) . ( 1731 and (1741) imply 

cos($ PR + j) = cos($ d + z/ cx vr + — + -) = (-1)"" cos($ d - -), (76) 

which explains the extra phase in (|68|) as well as the opposite signs on the 
7r/4 in the asymptotic formulas ([TBI and (j5_TJ) . 

In the classically forbidden region, the analytic continuations of the areas 
of the two lunes become complex, but their real parts are constant in any 
given region (ABCD), so only the imaginary parts need be equated. In this 
case, the condition is 

$pr = <W), (77) 

where $pr is given by (12011 and $<i by (|63l . Since both 3>pp and $^ vanish on 
the caustic, the condition ( 1771) implies that if the 6j-symbol is on a caustic, 
then so is the <i-matrix (as it should be). 
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To find the root of (1751) in the classically allowed region it helps to know 
the values of the areas of the lunes on the <i-sphere at the two turning points 
Pi and Pi, defined by flfJTj) . Call the half areas of the ci-lunes at these two 
turning points Ai and Ai- Then we have 

1 

A± = [j + — — max(m, m )]tt, A2 = max(0, — m — m )ir. (78) 

From these we compute an initial estimate for the root based on linear in- 
terpolation, 

fa = fa + Al -* PR f° (fa-fa), (79) 

after which a Newton-Raphson iteration converges to the actual root in all 
the cases we have examined, without taking iterations outside the interval 
[fa, fa)- We have no proof that this is always true, however. 

As for the classically forbidden regions, there are two of them, < (3 < fa 
and fa < P < ii, in which we must solve ( 177|) . Here we cannot use linear 
interpolation to find an initial estimate of the root since <3>d is not defined at 
P = or {3 = 7r. Instead we have simply taken 

n _ J \fa, < < fa, 

Po ~\ Kfa + ir), fa<P<ir, (80) 

and then applied a Newton-Raphson iteration. Occasionally this takes us 
outside the given classically forbidden region, whereupon we have reset the 
value of P to a point inside the region, using a simple prescription. After 
this, the Newton-Raphson converges to the root in all the cases we have 
examined. Our algorithm has proven satisfactory for our exploratory studies, 
but in more serious work the root finder will require more careful attention. 

The Newton-Raphson method requires us to know the derivatives d^d/dP 
or d^d/dp. These also appear in the role that $^ plays as a generating 
function, and are equal to the momentum pp of the rigid body, of which 
d 3 mm ,{P) is an eigenfunction (on the group manifold SU(2)). The derivatives 
are given by 

d$ d d$ d J\V d \ , 01 , 

or = . 81 

dp dp sin/3 V ; 

This applies at the intersection a (not a') in Fig. [T3J or at its analytic contin- 
uation as specified by the definitions of the angles ai at the end of Sec. 13.21 
It is convenient in using this formula to avoid calculating the sines of k, 
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or T], which become complex in the classically forbidden region; this can be 
done if |V^| is evaluated by taking the square root of the expression in (1551) . 

The value of (3 is smooth and well behaved as we cross from the classically 
allowed to forbidden regions, just as is the function X(x) in the method of 
comparison equations. This is illustrated in Fig. [TBI which uses the same 
parameters as Fig. [31 The heavy line in the figure is the same caustic line as 
in Fig. [31 and the light lines are contours of (3, labeled in degrees. 

Finally, we remark that the ratio of the amplitudes Ap^/A^ in (1681) has 
the form oo/oo as the caustic is approached, so a well designed numerical 
implementation of that formula would give special treatment to a small region 
around the caustic, where l'Hospital's rule would be used to avoid numerical 
difficulties. 

4.3 Numerical Results 

Figure 1T71 shows the results of numerical tests of the new uniform approxima- 
tion, with comparison with the Ponzano-Regge approximation. In the figure 
errors are plotted as a function of j'12 for the 6j-symbols 



in part (a) and part (b) of the figure, respectively. The values of the five fixed 
j's in part (b) are 8 times those in part (a). The plots show the absolute 
value of the difference between the exact 6j-symbol and the approximate 
value. In both parts of the figure, the curve labeled PR is the error of the 
Ponzano-Regge approximation, while that labeled U is the error of the uni- 
form approximation. The error of the Ponzano-Regge approximation is large 
near the caustics, as expected, while the error of the uniform approximation 
is fairly flat throughout the classically allowed region and up to the caus- 
tics. The error of both approximations falls rapidly to zero in the classically 
forbidden regions, as of course does the exact 6j-symbol. 

We computed the exact 6j-symbol with extended precision integer arith- 
metic. Recently other calculations have used extended precision floating 
point arithmetic to study the Wigner 3nj-symbols and their asymptotic 
properties [301 [31] . Extended precision is required when summing series with 
alternating (or variable) sign; often the desired sum is an exponentially small 
residue left when much larger terms nearly cancel. With floating point arith- 
metic one must carry enough extra precision so that enough remains after 




and 
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the subtractions; with integer arithmetic the answer is exact, but the can- 
cellations are still present (hence higher precision is carried on intermediate 
results). 

In the classically allowed region the error is oscillatory, and it is possible 
for the Ponzano-Regge error to be less than that of the uniform approxima- 
tion simply because it accidentally happens to fall near a zero of the cosine 
function. One such descending spike near j 12 = 140 can be seen in part (b) 
of the figure. It is clear, however, that a fair comparison of the errors in 
the classically allowed region must use the amplitude of the oscillatory func- 
tion and ignore the oscillations. By this measure the error of the uniform 
approximation in Fig. [T7| is approximately 30 times smaller than that of the 
Ponzano-Regge approximation in the center of the classically allowed region, 
and gets better as we approach the caustics. This ratio is nearly the same 
in parts (a) and (b), indicating that both errors scale in the same way with 
j. In the classically forbidden region the errors can be compared directly, 
without removing any oscillatory factor, and again for the values used the 
figure shows that the error in the uniform approximation is smaller than that 
in the Ponzano-Regge approximation. 

The error term for the Ponzano-Regge approximation is unknown, as 
is that for the uniform approximation, so there is no theory by which the 
errors can be compared. We would expect, however, on general grounds 
that the two error terms should scale the same with j, a conclusion that 
is supported by the numerical evidence. That the ratio between the errors 
should be as small as seen in fl82l) and Fig. [17] was a surprise to us, and we 
have no explanation for it. The values of the j's chosen in that example were 
essentially random, but when we try other "randomly chosen" values of the 
j's we get similar plots. If we systematically search for j values such that the 
error of the uniform approximation is as unfavorable as possible relative to 
that of the Ponzano-Regge approximation in the classically allowed region, 
we find cases such as 



which gives the error plots in Fig. [TS1 In this case the two errors are compara- 
ble for an extended range of j^. We have found no cases in which the uniform 
approximation is much worse than the Ponzano-Regge approximation in the 
classically allowed region. 

There is the question of when the uniform approximation is worst in an 
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absolute sense. We define the relative error as the difference between the 
approximation and the exact value, divided by a reference value. In the 
classically forbidden region, the reference value is the absolute value of the 
exact value. Inside the maximum part of the Airy function lobe around the 
turning points the reference value is the value of the Airy function, without 
any cosine modulation. Elsewhere in the classically allowed region we take 
the reference value to be the amplitude of the Ponzano-Regge approximation. 

Using this definition of relative error, we have systematically searched for 
j values that make the relative error in the uniform approximation largest. 
We find that they occur in cases for which ju = j'23 = 0. This can only 
happen when the other four j's are equal, so we have a 6j-symbol of the 
form 

f j j 1 

j 3 or- 



es*) 



For such symbols, the relative errors in both the uniform and the Ponzano- 
Regge approximation actually increase with j, reaching approximately 0.5 
(uniform) or 1.1 (Ponzano-Regge) when j = 10. The tetrahedra correspond- 
ing to 6j-symbols of the form (154"1) have two small edges (j'12 and j23; with 
lengths J12 = J23 = \) with no vertex in common, as illustrated in part (a) 
of Fig. M 

The uniform approximation is asymptotic, so it is no surprise that it does 
not work well for small quantum numbers such as those appearing in (l84lh 
Nevertheless it is interesting to see in more detail why the approximation is 
not good. Figure [201 shows the 6j-sphere for a symbol of the form (|84l . The 
orbit ji2 = (J12 = \) is the small circle about the south pole, while the 
orbit J23 = (J23 = |) is the curve that ends at the south pole in a cusp. 
The orbit J23 = is not smooth at the south pole, and it cannot be deformed 
into a small circle on the (i-sphere by any smooth map. Since the map is 
not smooth, it has infinite derivatives, and the semiclassical approximation 
breaks down. In fact, in this case, even the 6j-sphere itself is not smooth. 
That is, the 6j-sphere is obtained by symplectic reduction from a higher- 
dimensional phase space, and it is always topologically a sphere. But when 
the first four j's are all equal as in (1841) . the sphere is not differentiate at 
the south pole. In this case, it would be more appropriate to think of a tear 
drop with its cusp at the south pole, rather than a sphere. 

Unlike the Ponzano-Regge approximation, the uniform approximation is 
not invariant under all the symmetries of the 6j-symbol. It is invariant under 
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the three operations in which the upper and lower elements of two columns 
are swapped, but not under the six permutations of the columns. (We have 
not tested the "extra" symmetries due to Regge[2].) Therefore when j ^ 0, 
we can permute columns in an unfavorable case such as to obtain a better 
approximation. We have tested an algorithm in which, before applying the 
uniform approximation, the columns of the 6j-symbol are permuted to place 
the column with the largest minimum value in the third column. Then when 
we search for the worst case of the uniform approximation, we find that they 
occur with symbols of the form 



in which three of the j's are zero. Symbols of this form correspond to tetra- 
hedra that look like part (b) of Fig. [191 In this case the relative error in the 
uniform approximation no longer grows with j, rather it seems to approach 
a limit of about .075. 

Note that if any j in a 6j-symbol vanishes, the symbol can be evaluated 
trivially in closed form. Thus, no approximation is needed in the worst cases 
( |84|) and (|85|) that we have examined. 

5 Conclusions 

We have written computer codes that implement the uniform approximations 
given in this paper, which document the algorithms and check all their details. 
We will make these programs available to any interested parties. 

We will publish the derivation of the uniform approximation presented 
in this paper in a future article, in which we will make explicit the sym- 
plectic map between the 6j-sphere and the (i-sphere that underlies it, as 
well as outline how the theory of quantum normal forms leads to a uniform 
approximation in cases like this. 
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A Constructing the tetrahedron 

The 6j-symbol specifies the lengths Ji of the classical angular momentum 
vectors Jj but not their directions, so there is the question of how the ac- 
tual vectors can be constructed in three dimensional space, satisfying the 
identities (J2) and ©. 

Initially we assume that real vectors Jj exist, and we define 

Ai = Ji, 

A 2 = Ji+J2 = Jl2, 

A 3 = Ji + J 2 + J 3 = -J 4 , (86) 

which are the three vectors running along the edges emanating from the 
upper vertex in Fig. [2j We arrange these vectors as columns of a 3 x 3 
matrix F, and we let G = F T F, where T means transpose. Then G is the 
symmetric, nonnegative definite Gram matrix of dot products, 

Gij = A { ■ A r (87) 

By using the geometry of the three triangles spanned by the Aj, the compo- 
nents of G can be found in terms of the lengths Jj, 

Gii = A\ = J\ , G22 = A 2 
G\2 = G21 = Ai ■ A 2 : 

G13 = G31 = Al ■ A3 : 

G*23 = C32 = A 2 ■ A3 

Alternatively, without making any assumptions about the existence of the 
Ji, we can define G in terms of the lengths J, by (1551) . Then there is the 
question of whether vectors Aj exist such that (1571) is satisfied. 

The diagonalization of G is closely related to the singular value decom- 
position of F. The latter is F = UDV T , where U and V are real orthogonal 
matrices and D is a real diagonal matrix, containing the real singular values 
di on the diagonal. But this implies G = VD 2 V T , so V is the orthogonal 
matrix that diagonalizes G and the eigenvalues of G are df. Therefore to 
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2(^12 + J\ ~~ Ji) 
2(^1 + ^1 — ^23) 
5(^12 + J4 — J$) 



T 2 



(88) 
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find F we first construct G by (1881) and diagonalize it, obtaining V and the 
eigenvalues of G. If these eigenvalues are all nonnegative, then their square 
roots are the singular values, and the matrix D is determined. This does not 
determine U, but that matrix amounts to an overall rotation of the tetra- 
hedron which is arbitrary anyway. So we can set U to anything convenient, 
such as the identity. Then we have F = DV T , and the vectors Aj can be 
obtained as the columns of F . From these we can find the J's by inverting 
(J86D and using ©. 

Thus a real tetrahedron can be constructed if and only if the eigenvalues 
of G are nonnegative. The tetrahedron is only determined modulo overall 
rotations, proper and improper. If we wish the tetrahedron to have a definite 
handedness, we can perform a spatial inversion, if necessary, to make the 
volume in ffT3l) positive. The spatial inversion is properly brought about by 
time reversal, not parity, which does not change the sign of angular momenta. 
With this understanding, the final tetrahedron is determined modulo proper 
rotations. 

If any of the eigenvalues of G are negative, then a real tetrahedron does 
not exist. It turns out that at most one eigenvalue of G defined by ( 1551) 
can be negative, assuming the triangle inequalities on the Jj. A proof of 
this fact is easily given by considering the secular polynomial of G, whose 
coefficients can be expressed in terms of the lengths of the J's and the areas 
of the triangles formed by them. If G has one negative eigenvalue, then we 
can order the singular values so that the imaginary one is the third one. 
Then from F = DV T we see that the tetrahedron can be constructed with 
complex vectors, in which the x- and ^/-components of the J's are real, and 
the ^-components are purely imaginary. 
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Figure 1: The 6j-symbol as a function of ju for j\ = 16, ji = 80, jz = 208, 
j4 = 272 and j'23 = 276. Sticks are the values of the 67-symbol, and the curve 
is the Ponzano-Regge approximation. 

Figure 2: A tetrahedron of positive volume with conventional labeling of 
edges by angular momentum vectors. 

Figure 3: The J12-J23 plane for j 1 = §, j 2 = 3, j 3 = y, j 4 = 6. The 

Classical bounds are Jl2,min = f > ^12,max = \i ^23,min = |) <^23,max = y- ^ e 

dimension of the matrix (ji2 IJ23) is D — dimZ = 7. The point J 12 = 5, 
^23 = 9 (ji2 = § , J23 = y) is very close to the caustic line, but lies just 
inside. The Ponzano-Regge approximation is too large by a factor of 7 at 
this point. 

Figure 4: A sequence of four flat tetrahedra, moving around the caustic line 
of Fig. [3] in a clockwise direction from point Y. The parameters are the same 
as in Fig. [3j The numbers 1, 2 etc refer to vectors Ji, J2, etc. 

Figure 5: Definition of the interior dihedral angle 0i2- The other interior 
dihedral angle 023 is defined similarly. 

Figure 6: The phase space of the 6j-symbol is a sphere of radius D/2 in a 
space in which (K x , K y , K z ) are Cartesian coordinates. To within an additive 
constant, K z is J12 and the azimuthal angle 0i2 is the dihedral angle of the 
tetrahedron. Several curves of constant J 12 (small circles) are shown. 

Figure 7: Curves of constant J23 on the 6j-sphere. The first view shows the 
north pole and the point J 2 3 = ^23,min, and the second shows the south pole 
and the point J 2 3 = ^23,max- 

Figure 8: In part (a), the classically allowed region, an orbit of constant J12 
intersects an orbits of constant J 2 3- The shaded area is the Ponzano-Regge 
phase, to within an additive constant. In part (b), the classically forbidden 
region, the orbits do not intersect. 
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Figure 9: Caustics occur when the curve J\% = const is tangent to the curve 
J23 = const. The four types of such tangencies are illustrated. 

Figure 10: Definition of Euler angle (3 and unit vector n. 

Figure 11: Curves of constant J z and J n may intersect in the classically 
allowed region (a), or not intersect in the classically forbidden region (b). 

Figure 12: The square identifies the bounds on the classical observables J z 
and J n , while the spots indicate the quantized values J z = m, J n = ml '. The 
ellipse is the caustic curve. 

Figure 13: Vector a points to the intersection of the J 2 -orbit with the J n - 
orbit, with J y > 0. Vectors z, n and a define a spherical triangle, with 
interior angles 0, r\ and k. 

Figure 14: Caustics of the (i-matrices occur when the two small circles J z = 
const and J n = const are tangent. There are four possible configurations. 

Figure 15: In the method of comparison equations, the phase space of a non- 
linear oscillator (a) is mapped into the phase space of the harmonic oscillator 
(b). The function X(x) is determined by the equality of areas; for example, 
in the figure the shaded areas are equal, and X = X(xq). 

Figure 16: Contours of (3, the root of (!75|) or (ITT)) , in the J12-J23 plane. 
Parameters are same as in Fig. [3l but quantized values are omitted. 

Figure 17: Absolute value of the error of the Ponzano-Regge approximation 
(PR) and of the uniform approximation (U) as a function of j'12 for values 
of the other five j's shown in f)82p . The error is defined as the difference 
between the approximate value and the exact value. 

Figure 18: Comparison of Ponzano-Regge (PR) and uniform (U) errors as 
a function of j'12 for the 6j-symbol (|83|) . Dotted curve is the error in the 
uniform approximation. 
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Figure 19: The uniform approximation is worst for tetrahedra that look like 
part (a), where the small edges are ju and j'23. If columns are permuted to 
optimize the approximation, the worst case involves tetrahedra that look like 
part (b), where three edges are small. 



Figure 20: The orbits j'12 = is a small circle about the south pole, while the 
orbit J23 = has a cusp at the south pole. The shaded area is the "lune." 
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